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ABSTRACT 



Using the recently upgraded Long Baseline Array, we have measured the 
trigonometric parallax of PSR J0437-4715 to better than 1% precision, the most 
precise pulsar distance determination made to date. Comparing this VLBI dis- 
tance measurement to the kinematic distance obtained from pulsar timing, which 
is calculated from the pulsar's proper motion and apparent rate of change of or- 
bital period, gives a precise limit on the unmodeled relative acceleration between 
the Solar System and PSR J0437-4715, which can be used in a variety of ap- 
plications. Firstly, it shows that Newton's gravitational constant G is stable 
with time (G/G = (—5 ± 26) x 10~ 13 yr _1 , 95% confidence). Secondly, if a 
stochastic gravitational wave background existed at the currently quoted limit, 
this null result would fail ~ 50% of the time. Thirdly, it excludes Jupiter-mass 
planets within 226 AU of the Sun in 50% of the sky (95% confidence). Finally, 
the ~ 1% agreement of the parallax and orbital period derivative distances pro- 
vides a fundamental confirmation of the parallax distance method upon which 
all astronomical distances are based. 



Subject headings: astrometry — pulsars: individual J0437-4715) — gravitation 



1. Introduction 



PSR J0437-4715 has been observed intensively since its discovery by iJohnston et al. 
( 119931 ). It is t he brightest and nearest obser ved millisecond pulsar, and has also been studied 



in the optical (IBell et al.lll993l ). ultraviolet (IKargaltsev et al.ll2004l ). and X-ray (jZavlin et al. 
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20021 ) bands. The high rotational stability and close proximity of this pulsar-white dwarf 
binary system make it an excellent probe of General Relativity ( GR) and alternativ e form s 
of gravitational theories. The measurement of its Shapiro delay by Ivan Straten et al.l (120011 ). 
where the radio waves from the pulsar are delayed as they pass through the gravitational 
potential of the companion, is one such test which has shown consistency with GR predic- 
tions. The search for the low f requency stochast ic gravitational wave background (GWB) 
using pulsar timing arrays (e.g. IJenet et al.l 120051 ) is another test of GR which is facilitated 
in part by timing of J0437-4715. 

Although variation of Newton's gravitational const ant G with time is fo rb idden in GR, 
this is not required in alternate formalisms of gravity (e.g. lBrans fc Dickdll96ll ). IVerbiest et al 
(120081 ) have measured the rate of change of orbital period in the J0437-4715 system to 
better than 2% precision and shown that the agreement of the derived kinematic distance 
with the parallax distance derived from timing limits the time variation of G to less than 3 
parts in 10 11 . In this Letter, we present a new Very Long Baseline Interferometry (VLBI) 
determination of the position, proper motion, and parallax of PSR J0437-4715 and show 
that this improves the previously published G limit derived from this system by an order 
of magnitude, approaching the best publ ished limit (G/G = ( 4 ± 9) x 10~ 13 yr _1 ), which is 
derived from Lunar Laser Ranging (LLR; IWilliams et al.ll2004l ). Our VLBI observations and 
results are presented in §2] and ^respectively and the limitations on apparent accelerations 
(due to G or other possible causes such as unseen massive planets) are derived in £jH In ^5], 
we investigate the impact of the stochastic GWB on our results and derive an independent 
limit on the GWB amplitude. We summarize our conclusions in £JHJ 



VLBI observations 



Details of the observationa l, correlation, and post-processing strategies used are covered 
in detail in iDeller et al.l (120081 ) and only a brief summary is presented here. Four observa- 
tional epochs (MJD 53868, 54055, 54181, 54416) spanned 1.5 years, with each epoch lasting 
12 hours and utilizing all available antennas from the Australian Long Baseline Array (LB A). 
These were the three ATNF antennas (ATCA, Parkes, Mopra), the two University of Tas- 
mania antennas (Hobart and Ceduna) and, when available, one of the NASA Deep Space 
Network (DSN) antenna, at the Tidbinbilla statiotfl Observations were made at 8.4 GHz, 
and four 16 MHz bands were used, quantized at two bit precision to give a total data rate 
of 256 Mbps. 



1 DSS43 (70m) participated on MJD 54055 and 54181, while DSS34 (34m) was used on MJD 54416 
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Phase referencing was conducted using the International Celestial Reference Frame 
(ICRF) source J0439-4522, a bright (several hundred mJy flux at 8 GHz) and compact 



radio source less than 2° from PSR J0437-4715, whose position is known to 1 mas (IMa et al. 



19981 ). A six minute target/calibrator cycle was employed, with observing time divided 
equally between the target and calibrator. The well-studied blazar PKS 0537-441 was used 
for bandpass calibration. 

The data were correl ated, using matche d filtering (gating) on pulse profiles, with the 
DiFX software correlator (IDeller et al.l 120071 ). Due to the extremely narrow pulse profile of 
PSR J0437-4715 at 8.4 GHz, this gating improves the output signal to noise ratio (SNR) by 
a fac tor of 7 compared t o ungated data. The pulsar ephemeris used to set gates was taken 
from iHotan et all (bood ) . 



Due to the rel atively high frequen cy used for these observations, several of the calibration 
steps described in IDeller et al.l ( 120081 ) proved unnecessary. Data reduction was performed 
with and without GPS-based ionospheric correction, and the small (~ 100 //as) corrections 
which resulted from ionospheric correction gave a poorer parallax fit overall, with double 
the estimated systematic errors. Given the uncertainty of GPS-based corrections at these 
Southern latitudes, the ionospheric corrections were therefore not applied. Additionally, the 
effect of sci ntillation wa s negli gible at these higher frequencies, so the scintillation correction 
detailed in IDeller et al.l (120081 ) was not applied. 



Attaining single epoch accuracies of ~ 100 JJjdbS^ clS shown in §31 required post-correlation 
corrections to the visibilities due to the proper motion of the pulsar (~ 200 /ias in each 12 
hour observation). Additionally, the orbital motion of PSR J0437-4715 causes a peak to 
peak displacement of ~ 110 /ias. Compensation for this displacement in the single-epoch 
positions improved the quality of the proper motion and parallax fit, reducing the estimated 
systematic error contribution by 



/o. 



3. Astrometric results 



As discussed in IDeller et al.l ( 120081 ). the optimal weighting for visibility points when 
imaging the pulsar to fit for position depends on the ratio of thermal noise to residual sys- 
tematic errors. For strong targets systematic errors dominate the error budget, meaning 
the use of equally weighted visibilities (as opposed to weighting visibilities by baseline sen- 
sitivity) produces superior results. For PSR J0437-4715 the SNR of a typical single-band, 
single -epoch detection was ~20, well above the threshold of ~ 10 identified in IDeller et al. 
( 120081 ). and thus we would expect equally weighted visibilities to produce superior results 
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To confirm this, the parallax was determined from two separate datasets, which were pro- 
duced using sensitivity-weighted and equally weighted visibilities. As expected, the equally 
weighted dataset produced a better fit, with smaller errors on fitted parameters and smaller 
estimated systematic error contributions. 

The data were fitted using a n iterative approa ch designed to estimate and account 



for systematic errors, described in [Deller et al. 



shown in Table [T] along with the IVerbiest et al. 



arc 



2008). The resulting parameter values 
20081 ) timing measurement^; the VLBI fits 
to pulsar motion and observed positions are shown in Figure [TJ All errors are la unless 
otherwise stated, and include the covariances between fitted parameters. Over this short 
timespan, the covariance between parallax and proper motion is significant, and constitutes 
approximately half of the quoted parallax uncertainty. 

The average total single-epoch error of 132 //as is the best relative astrom etry performed 



by th e LBA to date, and is similar to other recent VLBI astrometric results (e.g lLoinard et al. 



20071 ). By way of comparison, the lim iting accuracy for EVN observations at 8.4 GHz has 
been simulated by iPradel et al.l (120061 ) to be 83 //as with a 1° calibrator-target separation, 
at a declination of 50°. The EVN has similar sensitivity to the LBA, but somewhat better 
uv coverage and longer baselines. 

Our derived distance of 156.3 ± 1.3 pc is the most accurate distance measurement (in 
both absolute and fractional distance) for a pulsar to date and approaches the most accurate 
distance measurem ents made of objects outside the solar system (T Tauri, 147.6 ± 0.6 pc: 
Loinard et al.ll2007l ). Previously, the high est-precision VLBI pulsar distance determinations 
were those made by iBrisken et al.l (120021 ). who measured the distance of PSR J0953+0755 
to an accuracy of 5 pc, along with 8 other Northern Hemisphere pulsars. Th e two previous 



paral lax measurements made using a Southern array, of PSR J0835-4510 (IDodson et al. 



2003J) and PSR J1456-6843 (IBailes et alJ Il990h . had la distance errors of 19 and 70 pc 
respectively. 



3.1. Comparison to timing astrometry 



To compa r e the VLBI and timing positions, we have re-fitted the timing data of 



Verbiest et al.l (120081 ) to obtain the position of PSR J0437-4715 at our reference epoch 



(MJD 54100). Table Q] shows that the timing and VLBI positions at MJD 54100 differ 



2 The timing positions of 
54100 



Verbiest et al 



(2008) have been re-fitted at the VLBI reference epoch of MJD 
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by over two mas, many times the formal errors shown. However, the formal errors are negli- 
gible compared to the 1 mas uncertainty in the VLBI phase reference calibrator position and 
the potential constant offsets due to phase-refer encing; errors such as station position errors 
(known to exist at the cm level for the LBA; iDeller et al.ll2008l ). Discrepancies between 
interferometri c and timing posit ions of even larger magnitudes have been found using the 
DE200 frame faartel et al.lll996h . and differences at the mas level still exist for the position 
PSR J0437-4715 calculated u sing the newer DE414 solar system ephemeris, as compared to 
the DE405 ephemeris used by lVerbiest et al.l (120081 ). Thus, we conclude that the VLBI and 
timing position difference is consistent with the uncertainty in the calibrator position and 
the offset between the solar system frame and the ICRF. 

The values obtained for the VLBI proper motion differ by ~ 4cr in both right ascension 
and declination from the timing values. A likely cause for this discrepancy is small changes 
in the centroid position of the phase reference source due to intrinsic source variability, 
which would be absorbed into the proper motion fit. The phase reference source, which is 
modeled by a ~ 1 mas FWHM Gaussian with additional delta components within 5 mas of 
strength 2% and 0.2% of peak flux, shows no gross evidence of variability (width/positions 
of primary /secondary components constant to ~ 200 /ias, and secondary fluxes are constant 
to ~ lmJy), but as it is only barely resolved (the beamsi ze is ~ 3 mas) variability at the 
~ 100 ^as level would be difficult to detect. iTitovl (120071 ) show that some ICRF sources 
exhibit apparent proper motions of hundreds of fias yr _1 for periods of several years, and 
while a detailed VLBI history of this source is not available, calibrator measurements from 
ATCA show that the flux density has declined by a factor of three over a five year period, 
implying some variability. The higher proper motion precision obtained with the timing data 
(a factor of 5-7 times better than the VLBI results) reflects the fact that the timing data 
spans a time baseline 7 times longer than the VLBI dataset. 

The parallax value obtained from VLBI is consistent with that derived from timing, and 
yields a distance which is consistent with the kinematic distance of 157.0 ± 2.4 pc derived 
from the orbital period derivative Pb- However, the VLBI parallax measurement is an order 
of magnitude more precise than the timing measurement, and yields a distance which is a 
factor of two more precise than the kinematic distance. 



4. Limits on anomalous accelerations 

The newly measured parallax of 7r = 6.396 ±0.054 mas allows an improved measurement 
of any anomalous acceleration of either the Solar System or PSR J0437-4715. Specifically, 
the apparent acceleration due to time variability of Newton's gravitational constant G and 
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the mass of an undetected trans-Neptunian planet near the line of sight to the pulsar can be 
limited. 



4.1. Constraints on G 



As first described by ID amour fc Taylor! (1199 ll ) , a precise measurement of a binary pul- 
sar's orbital period deriva tive, Ph, can be used to constrain a variation of the gravitational 



constant, G. However, as iBell &: Bailed (119961 ) pointed out, for PSR J0437— 4715 a precise 



distance needs to b e known in orde r to correct P^ for the Shklovskii acceleration caused by 
its proper motion ( S h klovskiil 1 1 9 701) . T his analysis has been performed, based exclusively 
on timing data, by 



Verbiest et al. 



(120081 ). whose limit was dominated by the uncertainty in 
their parallax measurement. Our VLBI parallax value, however, improves this limit by a 
factor of nearly 10 down to G/G = (—5 ± 26) x 10~ 13 yr _1 at 95% certainty. This value 
compares well to the most st ringent limit currentl y published: (4 ± 9) x 10~ 13 yr _1 , derived 



through lunar laser ranging (IWilliams et al.l 12004 ). Since 7r and P^ are now both measured 
to similar precision, both measurements will have to be improved for a further significant 
increase in G sensitivity. Additional VLBI observations and continued timing could see this 
limit improve upon the existing LLR limit early in the next decade. 



4.2. Acceleration due to massive bodies 



An alternative source of anomalous acceleration is he avy planets in a wide o r bit ar ound 
the Sun or the pulsar. Building upon the initial analysis of lZakamska fc Tremaind (120051). our 
parallax measurement can be combined with the timing results from I Verbiest et al.l (120081 ) 



to derive the follo wing result: a P Jc = (3 ± 16) x 10 s at the la level. This improves the 



limit published in I Verbiest et al.l (120081 ) by an order of magnitude and makes PSR J0437- 
4715 a more se nsitive Solar System accelerom eter than PSR J1713+0747, the most precise 
pulsar listed by Zakamska fc Tremaine ( 2005 ). From this, the limit for 50% of the skjj^l can 
be calculated as: |a©,5o%sk y /c| < 3.9 x 10 _19 s _1 (95% certainty). This acceleration limit 
can be used to exclude massive bodies within a given radius of the Sun; for example, at 
Kuiper-belt radii (50 AU) it excludes a planet more massive than Uranus over 50% of the 
sky, while Jupiter-mass planets are excluded within 226 AU over 50% of the sky. 



3 within 60° of the line of sight towards and away from PSR J0437-4715 
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5. Impact of the stochastic GWB 



Using tools recently developed by iHobbs et al.l (120081 ). we have simulated the effect 



of a GWB with spectral index —2/3 (as predicted for GWBs caused by black hole-black 
hole mergers) an d dimensionless amplitude 1.1 x 10~ 14 (the best published GWB limit; 



Jenet et al.l 120061 ) on the observed value of from pulsar timing. The simulated GWB 
causes the kinematic distance to be inconsistent with the VLBI parallax distance at the 
2cr level in ~ 50% of trials. Thus, although these observations cannot improve upon the 
present GWB limit, they are consistent with it. Simulations with a GWB with amplitude of 
1.1 x 10~ 13 show inconsistencies between the kinematic and VLBI distances at the 2a level in 
95% of trials, providing an independent exclusion of a GWB with an amplitude at or above 
this value. 

It is also interesting to note that the precise limit on G presented in §4.11 would be 
impossible in a Universe with a strong GWB. In the simulations with GWB amplitude 
1.1 x 10~ 13 , the observed G value is inconsistent with in 99% of cases, merely due to the 
GWB-induced corruption of the timing measurements. Thus, the stochastic GWB must 
eventually limit the accuracy of measurements of G in the fashion outlined in this paper. 



6. Conclusions 



We have observed PSR J0437-4715 using the LBA and obtained the most precise pulsar 
distance determination to date, with an error < 1.5 pc. Combined with accurate timing data, 
this has enabled us to confirm that \G/G\ < 3.1 x 10~ 12 yr -1 , the most stringent limit not 
obtained though Solar System tests. Alternatively, assuming an unchanging gravitational 
force, the results can be interpreted as excluding any unseen Jupiter-mass planets within 
226 AU of the Sun in 50% of the sky. Finally, the agreement between VLBI and timing 
results in this single case implies an upper limit to the stochastic GWB amplitude which 
is within an order of magnitude of the best limit derived from observations of ensembles of 
pulsars. 
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Table 1. Fitted VLBI results for PSR J0437-4715 and comparative timing values 
(positions re-referenced to the VLBI proper motion epoch) 



Parameter 


Fitted value and error 


Verbiest et al. (2008) ti 


ming values 


Right Ascension (J2000) 


04 h 37 m 15.883250 s ± 


0.000003 


04 h 37 m 15.883185 s 


± 0.000006 


Declination (J2000) 


-47° 15'09. 031863" ± 


0.000037 


-47° 15'09. 034033" 


± 0.000070 


Ma (mas/yr) 


121.679 ± 


0.052 


121.453 


± 0.010 


Mi (mas/yr) 


-71.820 ± 


0.086 


-71.457 


± 0.012 


Parallax 7r (mas) 


6.396 ± 


0.054 


6.65 


± 0.51 


Distance (pc) 


156.3 ± 


1.3 


157.0 


± 2.4 A 


Transverse velocity vt (km/s) 


104.71 ± 


0.95 


104.9 


± 1.6 A 


Reduced chi— squared 


1.0 








Average epoch mean fit error (mas) 


0.059 








Average intra— epoch systematic error (mas) 


0.068 








Average inter— epoch systematic error (mas) 


0.103 








Reference epoch for proper motion (M.TD) 


54100.0 




54100.0 





Derived from the kinematic distance obtained from , not the less precise parallax values 
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Fig. 1. — Motion of PSR J0437-4715, with measured positions overlaid on the best fit. 
Clockwise from top left: Motion in declination vs right ascension; as before but with proper 
motion subtracted; right ascension vs time with proper motion subtracted; and declination 
vs time with proper motion subtracted. 



